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The production of J/ty mesons in central collisions of heavy nuclei is investigated as a function of 
collision energy. Two contributions are considered simultaneously: early (hard) production coupled 
with subsequent suppression in a Quark-Gluon Plasma, as well as thermal recombination of primor- 
dially produced c and c quarks at the hadronization transition. Whereas the former still constitutes 
the major fraction of the observed J/ty abundance at SpS energies, the latter dominates the yield 
at RHIC. The resulting excitation function for the number of J/^'s over open charm pairs exhibits 
nontrivial structure around \fs ~ 30 AGeV, evolving into a significant rise towards maximal RHIC 
energy. We study this feature within different (thermal) scenarios for J /\P suppression, including 
parton-induced quasifree destruction as a novel mechanism. 



A promising probe for hot and dense QCD matter 
as created in the early phases of heavy-ion reactions is 
the abundance of J/\P mesons, measured via their de- 
cay branching into dilepton final states. Their produc- 
tion is expected to be significantly reduced in the case of 
Quark-Gluon Plasma (QGP) formation in sufficiently en- 
ergetic collisions of large nuclei. In such a picture, J/ty 
mesons are exclusively formed primordially (i.e., upon 
first impact of the colliding nucleons), and subsequently 
dissociated by (i) nuclear absorption, (ii) parton-induced 
destruction in a QGP |l| and/or Debye screening [^, and 
(iii) inelastic scattering on "comoving" hadrons in the fi- 
nal hadron gas phase of the reaction. The identification 
of the plasma effect thus requires a reliable knowledge of 
both (i) and (iii). Nuclear absorption is appreciable and 
has been thoroughly investigated in p-A and light-ion re- 
actions, whereas the impact of hadronic interactions is 
not yet well under control, although its net effect seems 
to be rather moderate, see, e.g., Refs. |3|-[7j . Within this 
framework the data of the NA50 collaboration at the SpS 
have been interpreted as evidence for QGP formation in 
the most central P&(158 AGeV)-P6 collisions [§. 

Recently, an alternative view of J/^ production in 
heavy-ion reactions has been put forward. Prompted by 
the observation that the J/*f? yield per charged hadron 
is remarkably constant as a function of impact param- 
eter, it has been argued in Ref. jjj] that all J/^P's are 
created statistically at the hadronization transition. The 
deduced temperature of T ~ 175 MeV is well in line with 
the so-called chemical freeze-out of light hadron produc- 
tion |]l0|-|l2||. However, no reference is made to an under- 
lying mechanism for cc creation. In a somewhat different 
approach, Braun-Munzinger and Stachel have extended 
their thermal model analysis Jio| to include (open and 
hidden) charm hadrons Jl3| (see Ref. jl4j for an update 
of this analysis). Together with the (dynamically well 
justified) proposition that cc pairs at SpS energies are 
exclusively produced primordially, the J/ty abundance 
in sufficiently central Pb-Pb collisions can be accounted 
for by statistical recombination of c and c quarks at the 
earlier inferred hadro-chemical freeze-out without intro- 
duction of new parameters. In a subsequent analysis [p~5|| , 



this approach was reiterated using a more complete set 
of charmed hadrons and enforcing exact (local) charm 
conservation within a canonical-ensemble treatment. Re- 
quiring to reproduce the NA50 measurements for J '/vp 
production, an open-charm enhancement factor of ~ 3 
relative to Np^ 3 times the value in N-N collisions was 
deduced for central Pb-Pb collisions (N p : number of par- 
ticipant nucleons). This coincides with the enhancement 
needed to explain the NA50 intermediate-mass region 
(IMR) dilepton spectra solely in terms of increased open- 
charm production^. As first pointed out in Ref. [ p^[ , 
the application of the thermal production framework to 
RHIC energies could in fact lead to an enhancement 
of J/^'s over its primordial production rate (see also 
Ref. H). 

One of the main assumptions in the thermal model 
analyses is that primordial production of J/'f's is absent, 
i.e., they either do not form or are fully suppressed before 
the hadronization transition. Under SpS conditions with 
supposedly rather short plasma lifetimes of 1-2 fm/c, this 
assertion is, however, not easily realized. In this article 
we therefore attempt a combined description of thermal 
and primordial production, the latter being subjected to 
nuclear absorption and plasma dissociation in an expand- 
ing fireball model around midrapidities. The focus will 
be on central Pb-Pb (Au-Au) collisions at varying CMS 
energy, covering the SpS and RHIC regime. We fix the 
participant number at N p ~ 360 (corresponding to an 
impact parameter b ~ 1.5 fm) to avoid complications as- 
sociated with transverse-energy fluctuations in the most 
central collisions |l^-pl[|. Also, cc pairs will be allowed 
to coalesce into charmonium states only over a limited 
range of rapidities. 

Let us start by discussing our implementation of the 
contribution from primordial production and subsequent 



*In this paper we adopt a scenario without any "anoma- 
lous" open-charm enhancement, attributing the NA50 IMR 
dilepton excess to thermal radiation Jl6|,[L7|. 
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suppression. As in Ref. (b|, we assume that the pro- 
duction of cc pairs entirely occurs in primary (hard) 
nucleon-nucleon collisions (secondary creation in a par- 
tonic medium has been shown to be negligible even at 
full RHIC energy p2[). In free space a certain fraction 
F<s, of these pairs combines into a final number NjL, of 
directly produced J/^f mesons (Fy ~ 2.5% in the SpS 
energy regime p3| ). In a heavy-ion environment, the 
first phase of suppression is characterized by interactions 
of the (pre-resonant) bound state with interpenetrating 
nucleons. It leads to a rather well understood (N p iN P 2) a 
suppression of the cross section with a = 0.92 ± 0.01 as 
inferred from p-A and A-B reactions with light projectile 
nuclei (see also Ref. p4| ). This factor equally applies to 
other charmonium states (\c, ^'), which contribute via 
their decay branchings into J/^f + X final states ("feed- 
down"). We model this so-called nuclear absorption us- 
ing a Glauber model with a phcnomcnological constant 
cross section <j^n ~ 5.8 mb pq| . 

In the second phase of suppression - the QGP - char- 
monium destruction has been discussed in both static 
screening-type pictures as well as dynamical ones via in- 
elastic collisions with partons, most notably the QCD 
analogue of photo-dissociation, g + J/ty — > cc |fj|j26f. 
Within an expanding fireball model (see below) we fol- 
low the dynamical picture, accounting, however, for a 
reduced J/\& binding energy. The dissociation rate is 
calculated from 
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with k m i n denoting the minimal on-shell momentum of a 
quark or gluon from the heat bath necessary to dissolve 
an in-medium charmonium bound state into a (free) cc 
pair. The binding is characterized by a temperature- 
dependent dissociation energy Edi SS taken from Ref. |27]] 
using a Debye screening mass m 2 D — g 2 T 2 with a typical 
g ~ 1.7; this entails E dlss (T = 180MeV) ~ 220 MeV, 
dropping to ~ 100 MeV at T = 240 MeV but crossing 
zero only around Toebye — 400 MeV. With such a de- 
crease in the J/\& binding energy, the break-up kinemat- 
ics render the photo-dissociation process increasingly in- 
efficient. For a loosely bound charmonium state, a more 
important mechanism turns out to be given by inelastic 
parton scattering, g(q, q) + J/^> — > g(q, q)+c+c. We eval- 
uate the respective cross sections in quasifree approxima- 
tion using leading-order QCD Q for gc — * gc (qc — > qc) 
and the appropriate break-up kinematics. In addition to 
a gluon screening mass, thermal quasiparticle masses for 
light quarks {ra 2 u d — g 2 T 2 /6, m 2 s — m\ + g 2 T 2 /6) and 
gluons (m 2 = g 2 T 2 /2) |29) are included. The resulting 

dissociation times, Tdi SS = r^L' are snown m Fig- an d 
compared to calculations with a constant break-up cross 
section of 1.5 mb, as well as to the photo-dissociation 
mechanism without medium effects in the J/ty bound 
state energy (as it has been employed in the literature 
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FIG. 1. Dissociation times for a J/ty in a QGP as a func- 
tion of temperature. The full curve corresponds to the lead- 
ing-order QCD process for quasifree g,q + c — > g,q + c scat- 
tering with in-medium J/ty bound state energy. The dashed 
curve represents a constant cross-section (Tdiss = 1-5 mb and 
the dotted curve results from photo-dissociation, gJ/^ — > cc, 
assuming the vacuum dissociation energy. 



before). At temperatures relevant for SpS conditions, 
the quasifree dissociation process (full line) is more ef- 
ficient than photo-dissociation (dotted line). At higher 
temperatures, it becomes less efficient due to an increas- 
ing gluon screening mass which suppresses the i-channel 
exchange graphs for g(q, q) + J/^ — > g(q, q) + c + c. 

The temperature (time) dependence of the dissociation 
rate has to be coupled with a model for the space-time 
evolution of the reaction dynamics. To facilitate the cal- 
culations we here employ a thermal fireball description 
in line with earlier analyses of dilepton radiation at both 
SpS (l6]fcj and RHIC |l)- Let us briefly recall its es- 
sential elements. After initial impact of two colliding Au 
(or Pb) nuclei the system is assumed to be thermalized 
after a formation time tq. Thereafter, the fireball un- 
dergoes isentropic expansion characterized by conserved 
entropy and (net) baryon number which defines a ther- 
modynamic trajectory in the [Ib — T plane of the phase 
diagram. Above the critical temperature T c a (quasipar- 
ticle) QGP equation of state is used, and a resonance 
hadron gas one below. The transition is modeled by a 
standard mixed phase construction |32| , 



S/V(t) = fs HG {T c ) + (1 - f)s QGP {T c ) , 



(2) 



justified for a sufficiently sharp increase of the entropy 
density around T c (/: fraction of matter in the hadronic 
phase, shg (sqgp) : entropy density in the hadronic 
(plasma) phase) . S denotes the total entropy in the con- 
sidered rapidity interval, and V(t) the time dependent 
volume therein, which we simulate by two fireballs with 
cylindrical expansion as 
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The parameters {w z ,a z ,a_L} are adjusted to finally ob- 
served flow velocities in connection with total fireball 
lifetimes of around 15 fm/c. The -y/s-dependence of the 
collisions is constructed as follows: the formation time To 
is taken to be ~ 1 fm/c at SpS and ~ | fm/c at full RHIC 
energy (with a powerlike interpolation in y/s), resulting 
in initial temperatures of T ~ 205 MeV and 390 MeV, 
respectively. We assume the transition temperature T c to 
smoothly increase from T c = 170 MeV at s/s — 17.3 GeV 
to T c = 180 MeV at ^/s = 200 GeV (with approxi- 
mately constant critical entropy density in the hadronic 
phase). SpS and first RHIC data (at ^/s = 56 GeV and 
130 GeV) on total multiplicities |33| and p/p, A/A ra- 
tios are used to estimate the total entropy as well 
as baryon and strange-quark chemical potentials. This 
description of the collision gives results consistent with 
hydrodynamical calculations, e.g., a pure QGP lifetime 
of ~ 1.5 fm/c (3.5 fm/c) at SpS (RHIC), and a mixed 
phase until ~ 5 fm/c (7.5 fm/c). The QGP suppres- 
sion factor Sqqp of J/ 1 ^ mesons follows from integrating 
the dissociation rate over the space-time evolution. The 
resulting direct yields per central collision (including nu- 
clear absorption) are listed in Tab. [jj 

The second source of charmonium states originates 
from thermal production at hadronization^. The under- 
lying picture Jl3|] is a statistical coalescence of c and c 
quarks at T c . In thermal models, hadron production is 
determined by the available phase space at T c . The total 
number of particle species j then is 
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where dj denotes the degeneracy factor, p,j the pertinent 
chemical potential, fij = Bjfis + SjfJ>s + Cj/i c , and V 
the hadronic fireball volume at T c . Since at CERN-SpS 
N cS <C 1, exact charm conservation is enforced within 
a canonical-ensemble treatment (see, e.g., Refs. f^J^l ) . 
We include all known charmed hadrons |37j] and fix the 
number N°£ r of cc pairs from primordial NN collisions 
in our restricted rapidity range as given by PYTHIA 
computations [p8[ upscaled by an empirical K factor, 



^In principle, J/ty formation can also occur above T c Jl^] 
through the reverse of the dissociation process, i.e., ccg — > 
J/^g; however, due to the smallness of the J/\E r binding en- 
ergy in the plasma, implying large formation times, we neglect 
formation above T c , thus possibly underestimating thermal 
production somewhat. 



K ~ 5 [|39|. This necessitates the introduction of a fu- 
gacity 7 C = j E for charm and anticharm quarks according 
to 
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where N open (Nhidden) denotes the thermal abundance of 
open (hidden) charm hadrons (Io,i are modified Bessel 
functions). Hence, the total contribution to statistical 
J/4" production (including strong and electromagnetic 
feeddown) follows as (J/*) = ^NfL, cf. Tab. | for se- 
lected collision energies. 
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TABLE I. Direct and statistical production of J/^'s per 
central collision (N p = 360) at various cm energies in a fixed 



rapidity window covered by two fireballs. AT 



number 



of primordial J/'f's remaining after nuclear absorption and 
plasma suppression (Sqgp). Jc- charm quark fugacity de- 
duced from ^ based on iV* r primordial cc pairs. (J/\f r ): 
number of statistically produced J/Ws . 

Combining the two sources of J/\I/'s (the suppressed 
direct as well as the statistical (or thermal) production), 
we calculate excitation functions from SpS to RHIC en- 
ergies. Fig. |^ displays the ratio of the observed number 
of J/^'s over the primordially produced one. The ther- 
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FIG. 2. Ratio of the number of observed J/ty's over 
the number of primordially produced ones (full curve) in 
a two-fireball model around midrapidity. The dashed 
(dash-dotted) curve corresponds to the J/^P yield from sta- 
tistical coalescence at hadronization (direct production with 
nuclear absorption and QGP suppression). 
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mal source accounts for one fourth of the yield at SpS, 
but dominates at higher energies implying a possible J/\E' 
enhancement at RHIC JT3| , p~8[ | . 

In anticipation of open-charm measurements at both 
RHIC and SpS E(J], one can make closer contact to ob- 
servables by plotting the ratio of the final number of 
J/^'s over the primordial number of cc pairs, cf. Fig. 
^. This ratio exhibits a minimum for y/s ~ 30 GeV when 
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FIG. 3. Ratio of the number of observed J/^'s over the 
number of primordial cc pairs (full curve) in a two-fireball 
model around midrapidity. The dashed (dash-dotted) curve 
shows the statistical contribution (direct production with nu- 
clear and QGP suppression) to this ratio. The band enclosed 
by the dotted lines reflects the uncertainty on the initial charm 
and J/ty production as explained in the text. 

the thermal and direct contributions are about equal. A 
similar minimum stucture has been found in Ref. pd|| , 
where, however, the J/\f r yield at all energies was entirely 
attributed to statistical production. 

The largest sensitivity in our calculation is attached 
to primordial cc production as indicated by the band en- 
closed by the dotted lines in Fig. [| The lower limit 
is estimated from next-to-leading order pQCD calcula- 
tions |4^] in connection with a lower bound in J/ty pro- 
duction from available data (supplemented by a phe- 
nomenological fit at higher energies ||43[| ). The upper 
limit is obtained from a PYTHIA calculation for cc pro- 
duction using GRV-HO structure functions (which tends 
to give the largest yield towards RHIC energies) in con- 
nection with J/vp production from Ref. Q. We also 
checked that there is only moderate sensitivity to vari- 
ations in the hadronization temperature: decreasing T c 
to 170 MeV at y/s = 200 GeV (with an accompanying 
increase in the hadronization volume, but at fixed N cS ) 
entails a 10% larger yield of thermal J/^'s (the decrease 
in the thermal density is overcompensated by the increase 
in volume and, more importantly, by the higher charm- 
quark fugacities). Within the uncertainties the plotted 



ratio persists to exhibit a very different behavior (i.e., an 
increase with CMS energy) from the one expected in the 
standard scenario of J/^> suppression. 

Finally, we investigate the sensitivity of the mini- 
mum structure with respect to different QGP suppression 
mechanisms. Upon replacing the quasifree destruction 
process by the gluon photo-dissociation process shown 
in Fig. |l| (dotted line), we observe a slight overall in- 
crease in the yield without significant alteration of the 
shape. Thirdly, in a more extreme scenario based on 
Debye-screening, we assume J/*f? mesons to be entirely 
suppressed if they are formed in a region with initial en- 
ergy density eo(r) > tDebye (along the lines of Ref. p4j). 
Within the Glauber model, the spatial distribution of 
primordial J/^'s is inferred from the nuclear thickness 
function Tab (?") (characterizing the number of A- TV col- 
lisions), whereas the energy-density profile is taken to be 
proportional to the density of participants in the trans- 
verse plane. We fix tDebye to obtain a suppression con- 
sistent with the NA50 data at y/s — 17.3 GeV (translat- 
ing into T Debye ~ 220 MeV). As expected, the pertinent 
excitation function exhibits a stronger suppression pat- 
tern with increasing ^/s, generating a more pronounced 
minimum structure (at similar position) in the ratios 
Nj /W /Nj% and N J/q ,/N^ than found with dynamical 
dissociation processes, cf. Fig. EL 
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FIG. 4. Comparison of Debye-screening (dashed curves), 
photo-dissociation (dotted curves) and quasifree dissoci- 
ation (full curves) of J/ty's in a thermalized QGP at 
various collision energies, combined with thermal produc- 
tion. Photo-dissociation (with vacuum J/$ binding en- 
ergy) gives results very similar to quasifree suppression (with 
in-medium Ed ias ) whereas Debye screening induces stronger 
suppression reflected in a more pronounced minimum around 
= 30 GeV. 

In summary, we have proposed a "combined" approach 
to evaluate J/^P yields in heavy-ion collisions which in- 
cludes (i) a direct contribution of prompt J/'J's sub- 
ject to nuclear and Quark-Gluon Plasma absorption and, 
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(ii) a thermal contribution of J/\&'s emerging from re- 
combination of c and c quarks at hadronization. The 
employed framework is consistent with earlier calcula- 
tions of thermal dilepton spectra; in particular, no en- 
hancement in open-charm production has been invoked. 
The resulting J/\I/ excitation function exhibits a transi- 
tion from mostly primordial to dominantly thermal pro- 
duction when going from SpS to RHIC. Such an inter- 
play could be mapped out by a systematic variation in 
collision energies accessible at RHIC. The predicted in- 
crease in J /ty yields will render J/^ suppression diffi- 
cult to identify as a QGP signature. However, the exci- 
tation function might serve as a sensitive probe of the 
hadronization dynamics at the QCD phase transition, 
provided an accurate knowledge of primordial cc abun- 
dances. We also note that our description might imply 
significant changes in the J/^f transverse-momentum dis- 
tributions. At the highest RHIC energies, one expects 
essentially thermal shapes (accompanied by a flow com- 
ponent from the QGP phase), which should be distin- 
guishable from hard production prevalent at the SpS. 

In this work, we did not address phenomena in the 
later hadronic stages of the collision, e.g., possible con- 
sequences of (the approach towards) chiral symmetry 
restoration via in-medium modifications of Z3-meson 
masses which might play an important role for the 
v^'/vj/ ratio |45| , [46| |. An extension of our approach along 
these lines together with a detailed comparison to avail- 
able data on centrality and projectile dependence is in 
progress [E7J. 
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